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Here we show the connection between topological order and the geometric entanglement, as mea- 
sured by the logarithm of the overlap between a given state and its closest product state of blocks, 
for the topological universality class of the toric code model. As happens for the entanglement 
entropy, we find that for large block sizes the geometric entanglement is, up to possible subleading 
corrections, the sum of two contributions: a non-universal bulk contribution obeying a boundary 
law times the number of blocks, and a universal contribution quantifying the underlying pattern of 
long-range entanglement of a topologically-ordered state. 



Introduction.- Topological order (TO) [1] is an exam- 
ple of new physics beyond Landau's symmetry-breaking 
paradigm of phase transitions. Systems exhibiting this 
new kind of order are linked to concepts of the deepest 
physical interest, e.g. quasiparticle anyonic statistics and 
topological quantum computation [2], Importantly, TO 
finds a realization in terms of Topological Quantum Field 
Theories [ ] , which are the low-energy limit of quantum 
lattice models such as the toric code and quantum double 
models [4] , as well as string-net models [5] . 

A remarkable property about TO is that it influences 
the long-range entanglement in the wave function of the 
system. For instance, as proven for systems in two spa- 
tial dimensions (2D) [6, 7], the entanglement entropy 
S of a block of boundary size L ^> 1 obeys the law 
S = So — 5 7 + 0(1/L), where Sq oc L is a non- universal 
term (the so-called "boundary law"), and S*y is a uni- 
versal long-distance contribution: the topological entan- 
glement entropy. 5 7 is non-zero for systems with TO, 
e.g. Sy = 1 for systems in the topological universality 
class of the toric code. More recently, similar univer- 
sal contributions have also been found for other bipartite 
entanglement measures [8, 9]. 

The purpose of this letter is to investigate, in the con- 
text of systems with TO, a global measure of entangle- 
ment which captures truly non-bipartite correlations in 
the system: the geometric entanglement ( GE) [1 0], which 
we call Eq. For simplicity, we focus on systems in the 
topological universality class of the toric code [4]. Our 
results indicate that, generically, the geometric entangle- 
ment of blocks of boundary size L>1 obeys the law 

E G = E -E y + O{l/L), (1) 

where E cx n^L is a non-universal term that obeys some 
boundary law scaling multiplied by the number of blocks 
rib, and E~ 1 is the topological geometric entanglement. 
Moreover, we find that E 1 = Sj, with S* 7 the topologi- 
cal entanglement entropy. These results are also the first 
explicit example of boundary laws for the GE in ground 
states of 2D quantum many-body systems. 



The toric code topological universality class.- Let us 
start by reviewing some basics on the toric code model 
[ ]. This is the renormalization group (RG) fixed point 
of the topological universality class of a Z2 lattice gauge 
theory, and is equivalent under local transformations to 
the Levin- Wen string model on a honeycomb lattice [5] 
(see the Appendix). 

We consider a square lattice E on a torus. Other Rie- 
mann surfaces of genus g could also be considered eas- 
ily without changing our conclusions. Non-bipartite lat- 
tices (e.g. honeycomb) could also be considered without 
changing the long-distance properties. There are spin- 
1 /2 (qubits) degrees of freedom attached to each link in 
E. The model is described in terms of stars and pla- 
quettes. A star "+" is a set of links sharing a common 
vertex. A plaquette "□" is an elementary face on the lat- 
tice E. For any star + and plaquette □, we consider the 
star operators A + and plaquette operators Brj defined 

as A + = rije+ cr ^ and b O = EL/en "*') where CTq 1 is 
the a-th Pauli matrix at link j of the lattice. Let us call 
respectively n + , n\j and n the number of stars, plaque- 
ttes and links in lattice E. Star and plaquette operators 
satisfy the global constraint 

]jA + =l[B a =I. (2) 

+ □ 

Therefore, there are n+ — 1 independent star operators 
and n\j — 1 independent plaquette operators. 

With the definitions above, the Hamiltonian of the 
model reads H = — J2 + A + — J2a^a- This Hamilto- 
nian can be diagonalized exactly as explained in Ref. [4] . 
The ground level of H is 4-fold degenerate (4g-fold for a 
Riemann surface of genus g). This degeneracy depends 
on the underlying topology of E, and is a signature of 
topological order. Moreover, the ground level is a stabi- 
lized space of G, the group of all the possible products of 
independent star operators, of size \G\ = 2(™+~ 1 ). 

In order to build a basis for the ground level subspace, 
let us consider a closed curve 7 running on the links of 
the dual lattice E*. We define its associated loop oper- 
ator W x py] = rijG 7 <J ^, where j £ 7 are the links in E 
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crossed by the curve 7 connecting the centres of the pla- 
quettes. It is not difficult to see that the group G can 
also be understood as the group generated by all con- 
tractible loop operators [23]. Let also 71 and 72 be the 
two non-contractiblc loops on a torus, and define the as- 
sociated string operators = W x [ji^]- We call |0) 
and |1) the eigenstates of a z respectively with +1 and 
— 1 eigenvalue. The ground level subspace then reads 
C = spanjji, j), i,j = 0,1}, where 

M = Kr 1/2 £s<^ ior n . (3) 

g£G 

It is easy to check that the four vectors \i,j) are orthonor- 
mal and stabilized by G, that is, g\i,j) — \i,j) Vg e G 
and Vi, j. These four states form a possible basis of the 
ground level subspace of the model on a torus. 

Excited states can be constructed by locally apply- 
ing Pauli operators a x , a v , a z on the ground states 
Pauli operators a z create pairs of deconfined charge- 
anticharge quasiparticle excitations, o~ x create pairs of 
deconfined flux-antiflux quasiparticles, and o~ y creates a 
flux-antiflux and a charge-anticharge pairs. These oper- 
ators can be applied to several sites of the lattice, thus 
creating a quasiparticle pattern that defines the excited 
state. The excited states of the model are then labeled by 
a set of quantum numbers, \<f>, c, where 4> and c are 
patterns denoting the position of flux-type and charge- 
type excitations respectively, and i,j label the ground 
state that was excited. In this paper we will focus 
on the entanglement properties of states \<j>, c, 

A key property of the toric code model is that its 
ground state |0,0) = \i = Q,j = 0) can be created by 
a quantum circuit that applies a sequence of Controlled- 
NOT (CNOT) unitary operations over an initial separa- 
ble state of all the qubits [11]. This means that it is actu- 
ally possible to disentangle qubits from the ground state 
of the system by reversing the action of these CNOTs. 
This was the key observation that allowed to build an 
exact MERA representation of the ground states \i,j) 
of the model [12]. The two fundamental disentangling 
movements are represented in Fig l(a,b), and leave the 
overall quantum state as a product state of the disentan- 
gled qubits with the rest of the system. What is more, 
the rest of the system is left in thejjround state of a toric 
code model on a deformed lattice S, where £ is obtained 
from E by removing the links that correspond to the dis- 
entangled qubits. This property turns out to be of great 
importance for some of the derivations in this paper. 

Geometric entanglement.- Let us now remind the ba- 
sics of the geometric entanglement. Consider an Tri- 
partite normalized pure state € % = ®fcLi% j 
where W-^ is the Hilbert space of party i. For instance, in 
a system of n spins each party could be a single spin, so 
that m = n, but could also be a set of spins, either con- 
tiguous (a block [13]) or not. We wish now to determine 




FIG. 1: (color online) (a,b) CNOT operations that disentan- 
gle qubits from the system, thus removing their links from 
the lattice. In the diagram, the arrows go from controlling to 
target qubits; (c) Example of a partition into blocks of a given 
boundary L, e.g. L — 12 here (but all our derivations work for 
arbitrary L). Qubits crossed by the boundary of a block are 
regarded as being inside the block; (d,e) Doing CNOT oper- 
ations locally inside of each block we can remove all the stars 
inside of all the blocks in two steps: first, we apply CNOTs 
as in (a) to get (d) and, second, we apply CNOTs as in (b) 
to get (e). Notice that in (e) each star is made of either 4 or 
8 qubits. The procedure works similar for other blockings. 

how well state Y$>) can be approximated by an unentan- 
gled (normalized) state of the parties, |$) = <&™ =1 | 
The proximity of \^f) to |$) is captured by their overlap. 
The entanglement of |W) is thus revealed by the maximal 
overlap [10], A max (*) = max* |($|*)|. The larger A max 
is, the less entangled is \*S?). We quantify the entangle- 
ment of |\&) via the quantity: 

£ G (*) = -log 2 AL x W, (4) 

where we have taken the base-2 logarithm, and which 
gives zero for unentangled states. Ec{^) is called geo- 
metric entanglement (GE). This quantity has been stud- 
ied in a variety of contexts, including critical systems 
and quantum phase transitions [13, 14], quantification 
of entanglement as a resource for quantum computation 
[15], local state discrimination [16], and has been recently 
measured in NMR experiments [17]. Also, one can choose 
the case of just two sets of spins. In this case the GE 
Ec(^) coincides with the so-called single-copy entangle- 
ment between the two sets, £'i(4 r ) = — log 2 fi(p), with 
v\ (/?) the largest eigenvalue of the reduced density matrix 
p of either set [18]. 

The GE offers a lot of flexibility to study multipartite 
quantum correlations in spin systems. For instance, one 
can choose each party to be a single spin, but one can also 
choose blocks of increasing boundary length L [13] (see 
Fig. 1(c) for an example). Studying how the GE changes 
with L provides information about how close the system 
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FIG. 2: (color online) Bipartite lattices of spins and blocks: 

(a) Sets of spins A and B for Theorem 1 and a 4 x 4 torus; 

(b) Sets of blocks A and B for Theorem 2. 

is to a product state under coarse-graining transforma- 
tions. 

GE of the toric code universality class.- We now study 
the GE of the topological universalty class of the toric 
code model. By definition of topological universality 
class, the long-distance properties of all the states in this 
class are equivalent to those of the RG fixed point, that 
is, the toric code model. Hence, if we are interested in ex- 
tracting the universal topological (long-distance) contri- 
bution to the GE of the states in this class, it is sufficient 
to focus our analysis on the fixed point only. 

We thus consider the entanglement properties of the 
four ground level states \i,j) as well as the excited states 
\4>, c, Our main result here is that, for these states, 
the GE of blocks consists of a boundary term plus a topo- 
logical contribution. In order to see this we first consider 
a couple of useful Lemmas: 

Lemma 1: All the eigenstates \<j>, c, of the toric 
code Hamiltonian have the same entanglement properties. 

Proof: let us choose a reference state from the ground 
level basis, e.g. |0, 0). First of all, notice that the four 
basis states \i,j) in the ground level are related to |0, 0) 
by local unitary operators, since the string operators w\ 
and W2 are tensor products of 2 x 2 identity operators 
and Pauli-x matrices. Since these operators act locally 
on each spin, they do not change the entanglement of the 
quantum state. Moreover, the excited states \<f>, c, 
are created locally by applying tensor products of 2 x 2 
identity and Pauli-x and z operators to \i,j). The states 
resulting from these operations have then the entangle- 
ment properties of state \i,j), which in turn are the same 
as those of the reference state |0,0). Thus, all states 
\4>, c, have the same entanglement properties. □ 

Lemma 2: Consider an arbitrary set of blocks of 
qubits in lattice S, where each block is regarded as an in- 
dividual party. Then, the entanglement properties of the 
ground state |0, 0) are the same as those of state |0,0), 
the ground state of the toric code model in the deformed 
lattice S obtained after disentangling as many qubits as 
possible using CNOTs locally inside of each block. 

Proof: given the ground state |0,0) in S, we start 
by doing CNOT disentangling operations locally inside 



of each block in order to disentangle as many qubits as 
possible. As a result, state |0,0) transforms into state 
1 0, 0) disentangled = |ei) ® ••• ® \e p ) ® |0,0), where \e k ) 
is the quantum state for the k-th disentangled qubit, 
k = l,...,p (p is the number of disentangled qubits), 
and 1 0, 0) is the i = 0,j = ground state of a toric 
code Hamiltonian in the deformed lattice S. Since all 
CNOT operations are done locally inside of each block, 
states |0,0) and |0, 0)di S entangicd have the same entangle- 
ment content if we regard the blocks as individual parties. 
What is more, the entanglement in the system is entirely 
due to the qubits in |0, 0), which proves the lemma. □ 

Keeping the above two lemmas in mind, we now 
present the following theorem about the geometric en- 
tanglement of spins in the toric code: 

Theorem 1 (for spins): For the toric code Hamilto- 
nian in lattice S, the four ground states for i,j — 
0, 1 and also the excited states \(f>, c, i,j) have all the same 
GE of spins and is given by 

E G = n+ - 1, (5) 

where n + is the number of stars in S . 

Proof: Using Lemma 1 we can restrict our attention to 
the GE for the state |0, 0). In this setting, we call compu- 
tational basis the basis of the many-body Hilbert space 
constructed from the tensor products of the {|0), |1)} lo- 
cal basis for every spin, which is an example of product 
basis for the spins. 

Our proof follows from upper and lower bounding the 
quantity E G (0,0) in Eq.(4) for the ground state |0, 0). 
First, from the expression in Eq.(3) for the ground states 
we immediately have that the absolute value of the over- 
lap with any state of the computational basis is |G|~ 1//2 , 
and therefore A max > | C| x / 2 . From here, we get 
E G (0,0) < log 2 \G\, which gives an upper bound. 

Next, to derive a lower bound we use the fact that 
if we group the n spins into two sets A and B, then 
A max < Am a x ' , where [A : B] means that a partition 
of the system with respect to the two sets A and B is 
considered. Thus, we have that E G (0,0) > E G (0,0)^ A:B l 

The trick to find a useful lower bound is to find an ap- 
propriate choice of sets A and B. In our case, we consider 
e.g. the bipartiton shown in Fig. 2(a) for even x even lat- 
tices (other cases can be considered similarly) . Then, we 
use a Lemma by Hamma, Ionicioiu and Zanardi in Ref . [ ] , 
that the reduced density matrix pa of |0, 0) for subsys- 
tem A satisfies p\ — (\G a\\Gb\/\G\)pa, where Ga/b is 
the subgroup of G acting trivially on subsystem A/B. 
Importantly, for A and B chosen as explained above, it 
happens that G a and Gb are trivial groups consisting of 
only the identity element. With this in mind, we see that 
the reduced density matrix pa has eigenvalues either zero 
or |G| _1 , which is |G|-fold degenerate. This means that 
(A^f 1 ) 2 = IG]" 1 and hence E G (0,0) > log 2 |G|. 
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Combining the two bounds, we get Eg (0,0) = 
log 2 \G\ = n + — 1, and from here Eq.(5) for the GE for 
spins follows immediately. □ 

Importantly, the techniques used in Theorem 1 can also 
be used to deal with blocks of spins whenever the blocks 
form a bipartite lattice. In general, one may first disen- 
tangle as many qubits as possible inside the blocks. Then 
the remaining qubits are left in an entangled state where 
the GE is equal to the number of remaining independent 
star operators, which amounts to a boundary law term 
(possibly with a subleading correction depending on how 
the blocks are chosen) plus a topological term. As an 
example of this, let us present the following theorem: 

Theorem 2 (for blocks): Given a partition of the 
lattice E into rif, blocks of boundary size L as indicated 
in Fig. 1(c) (where L is measured in number of qubits), 
then the four ground states \i,j) for i,j — 0, 1 and also 
the excited states \<f>, c, of the toric code Hamiltonian 
have all the same GE of blocks and is given by 

E G = (n b L/i) -I. (6) 

Proof: First, notice that Lemma 1 and Lemma 2 imply 
that we can entirely focus on the ground state |0, 0) of a 
toric code model in the deformed lattice E from Fig. 1(e). 
As shown in Fig. l(d,e), it is always possible to remove all 
the stars inside of each block in E just by doing CNOTs 
locally inside of each block. As a result, the stars in 
the deformed lattice E correspond to those in E that lay 
among the blocks. It is easy to see that, for nj, blocks of 
boundary L, there are n+ = n^L/A of such stars. 

The rest of the proof follows as in Theorem 1, by up- 
per and lower bounding Eg(0, 0) for the partition with 
respect to blocks. For the upper bound, we use the fact 
that A£i° x cks] > A& ] > iq- 1 / 2 , where the first in- 
equality again used the property that if larger blocks are 
considered then the maximum overlap is also larger, and 
G is the^corresponding group of contractible loop opera- 
tors on E. From here Eg(0, 0) < log 2 \G\ follows. For the 
lower bound, we divide the system into two sets A and 
B of blocks as indicated in Fig. 2(b). Following a simi- 
lar reasoning as in Theorem 1, it is possible to see again 
that no element g £ G will act trivially on A or B except 
for the identity element. From this point, the rest of the 
proof is simply equivalent to the proof for Theorem 1. □ 

Discussion.- Let us now discuss Eqs. (5) and (6) in 
detail. First, let us remind that a variety of works have 
shown the existence of a "boundary law" for the entangle- 
ment entropy of many 2D systems [19], including models 
with TO [6, 7]. It is then remarkable that, according 
to the first term of these equations, the entanglement 
per block obeys also a boundary law. To the best of 
our knowledge, these results are the first example of a 
boundary law behaviour for a multipartite (rather than 
bipartite) measure of entanglement in 2D. 



The second term in Eqs. (5) and (6) is far more intrigu- 
ing and important. Its existence is caused by the global 
constraint from Eq. (2) on star operators which, in turn, 
allow for the topological degeneracy of the ground state. 
Thus, this term is of topological nature and universal, and 
quantifies the pattern of long-range entanglement present 
in topologically ordered states. 

In order to clarify further the meaning of the topo- 
logical term, let us consider the case of just two blocks 
of spins. In this case, the geometric entanglement 
Eg coincides with the single-copy entanglement E\ = 
— log v\ (p) , with ui (p) the largest eigenvalue of the re- 
duced density matrix p of either block. Since this density 
matrix is proportional to a projector (see Ref. [ ] and also 
Ref. [9]), we have that E\ = S, with S the entanglement 
entropy of the bipartition. Thus Eg = E\ = S. We also 
know that for a system with TO, the entanglement en- 
tropy satisfies S — So — S 7 , where So is some boundary 
law term, and S" 7 is the topological entropy. As explained 
in Ref. [6] , for the toric code model the global constraints 
in Eq.(2) imply that 5* 7 = 1. Thus, with the conven- 
tion from Eq.(l), we have that the topological geometric 
entanglement is given by E 7 = 5 7 . 

Robustness of _E 7 away from the fixed point.- There 
are several ways of checking the robustness of E-y under 
perturbations. One possibility is to perform a numerical 
analysis. This is not considered here, and will be the 
subject of a future work [ ]. Another option is a pertur- 
bation theory analysis, which we present in detail in the 
Appendix. Yet, a more intuitive alternative is the follow- 
ing argumentation based on RG fixed points: the toric 
code on a square lattice is an RG fixed point and hence 
representative of its topological universality class. The 
long-distance properties of any state in this universality 
class do not change under local RG transformations and, 
thus, are equivalent to those of the fixed point. This, in 
particular, is true for the long-range pattern of entan- 
glement, and hence for E 1 . Thus, any non-relevant and 
short-range perturbation driving the Hamiltonian away 
from the fixed point will produce ground states with the 
same E~. Nevertheless, we expect a change in the short- 
range pattern of entanglement, and hence in the non- 
universal boundary law for Eq. For short-range pertur- 
bations the change involves modifications of the prefactor 
of the boundary law as well as the possible appearance 
of subleading 0(1/ L) corrections. This can be checked 
e.g. by perturbation theory (see the Appendix). Hence, 
Eq.(l) applies when away from the fixed point. 

Conclusions.- Here we have studied the GE of blocks 
in topologically-ordered states. We found that it is com- 
posed of a boundary law term plus a topological term, 
apart from possible subleading corrections. We focused 
on the topological universality class of the toric code 
model, which includes the simplest string-net model of 
Levin and Wen [ ] . We believe that similar results should 
also apply to the A-phase of Kitaev's honeycomb model 
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FIG. 3: (color online) Mapping of the toric code from a square lattice E to a honeycomb lattice E (Levin- Wen string model), 
and back to a square lattice E'. Each red arrow represent a CNOT operation that disentangles qubits from the system, as 
explained in Fig.l(a,b) of the main paper. The arrows go from controlling to target qubits. 



[21] for which the Toric Code is (in some limiting cases) 
an effective model, as well as spin-liquid states with an 
emergent Z 2 gauge symmetry [22]. 
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Appendix 

In this Appendix we provide the following further in- 
formation: first, we explain an explicit local mapping 
that shows the local equivalence of the toric code model 
on a square lattice and the Levin- Wen string model on 
the honeycomb lattice. Hence, bot models are valid RG 
fixed points representing the same topological univer- 
sality class. Second, we perform a perturbation theory 
analysis of the robustness of the GE of the toric code 
model under external short-range perturbations such as 
magnetic fields. This analysis provides upper and lower 
bounds for the GE, and complements the RG argumen- 
tation on the robustness of E-y presented in the paper. 



(i) Local mapping between the square and 
honeycomb lattices 



fixed point, see e.g. the Entanglement Renormalization 
analysis from Ref.([12]). 



(ii) Perturbation theory analysis of the robustness 

Let us now add a perturbation to the toric code Hamil- 
tonian on the square lattice, and see how the ground state 
changes. For simplicity, we consider the case of an infinite 
plane, where the ground state 1 , 0) is non-degenerate, 
and hence we can use non-degenerate perturbation the- 
ory. The perturbed Hamiltonian will be 



H x = H + XV, 



(7) 



where H is the toric code Hamiltonian, V is the pertur- 
bation, and A<1. Non-degenerate perturbation theory 
says that the new ground state can be approximated as 



|0,0) A «|0,0) 



^ (&c, 0,01^10,0) ,^ „ ... 
X l^— B 16 c >°>°) 



0,0 



E. 



<f>,c 



(8) 



where E is the ground state energy and E^ c is the 
energy of the excited state \<j>, c, 0, 0). 

We now consider the case in which the perturbation is 
an homogeneous magnetic field, e.g. in the x direction, 



The toric code model on a honeyconb lattice is, in fact, 
the simplest of the string-net models of Levin and Wen 
[5]. By definition, this model describes the fixed point of 
a topological phase capturing all the long-range proper- 
ties of the universality class of a Z 2 lattice gauge theory. 
Here we show that, in fact, the toric code model on a 
square lattice is also a valid RG fixed point. A way to 
see this is to realize that it is actually possible to map the 
model on the square lattice to the honeycomb lattice, and 
back again to the square, by means of a sequence of disen- 
tangling (coarse-graining) CNOT transformations acting 
locally on the lattice. These are represented in the dia- 
gram of Fig. (3). We stress, though, that it is already well 
known that the toric code in a square lattice is an RG 



f = 5>£' ] ( Q ) 

(the case of z and y directions can be considered simi- 
larly). It is easy to check that in this case, the normalized 
perturbed ground state becomes 

|o,o) A «c||o,o)- ^E^ ] l°'°)j ' ( 10 ) 

with A the energy gap to create a pair of flux and anti- 
flux quasiparticles, and C = (1 + nA 2 /A 2 )~ 1 / 2 a normal- 
ization constant. 
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Our aim now is to estimate the maximum overlap 
of the pervious state with a product state of blocks of 
boundary size L. This can be done as follows: first, and 
as in the unperturbed case, we apply CNOT operations 



locally inside of the blocks so that qubits are disentan- 
gled in the unperturbed ground state. By doing this, we 
can focus on the entanglement of the state 



10, 0)^ iscntanglcd « C I |0, 0} - ±fl a U] |o, 0) I ® | ei ) ® • • • ® |e p ) , (11) 



A 



where je^} is the quantum state for the fc-th disentangled qubit. The above equation is indeed equivalent to 



A " b ~ A ~ 

|0,0)^ cntanglcd «C ( |0,0> ® let) (8, - - . (g, | e?J ) - — ^^10,0) O | ei > o - - . ® |e p > - ^-|0, 0>lo; 1; ... !P ) ] , (12) 

I 



where Si? is the total spin in the a; direction for the L 
spins in the boundary of block j, and 

p 

K..., P >=^crM|e 1 )®---®|e p ) . (13) 

Now we find upper and lower bounds to the maximum 
overlap of this state with a product state of the blocks. A 
lower bound can be easily obtained by the product state 
|0}®(™-f) ® | ei ) ® • • • ® |e p ). Noticing that |e fe ) is either 
|0) or |+) [12], and that the |+) contributions come only 
from qubits close to the boundary of the block, we have 
that 

C^-^(l-^)<AK^, (14) 

for some positive u> = 0(1) constant. The following up- 
per bound can also be found easily: 

AN? x cksl < C\G\^ (l + ^ + A) . (is) 

Using the above bounds, one can check that for the GE 
we obtain, in the limit A <§; A and L ^> 1, 

(\ 2uj\\ , (\ 2A\ , . 

(4 + "A" j " fcL " 1 " ^ " (4 - A j - 1 ' (16) 

The above equation is compatible with a leading change 
in the GE in the prefactor of the boundary law. Also, 
the fact that both bounds leave the topological compo- 
nent Ery = 1 untouched seems to indicate that this is 
actually robust under the perturbation. Moreover, im- 
plementing finite- L corrections to these bounds provides 
0(1/ L) corrections, and thus Eq.(l) in the main paper 
applies. 
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